Multi-photon decays of the Higgs boson in standard model: 
leading terms from Heisenberg-Euler effective Lagrangian 
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Abstract 

We calculate the multi-photon decay widths of the Higgs boson from an effective La- 
grangian for a system of electromagnetic and Higgs fields. We utilize a low-energy theorem to 
connect the above Lagrangian to the Heisenberg-Euler effective Lagrangian induced by charged 
particles that gain mass from interactions with the Higgs boson. In the standard model only the 
W ± gauge bosons and the top quark are relevant, and we compute their contributions to the 
effective couplings and the multi-photon decay widths of the Higgs boson. 
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A new boson has recently been discovered at the CERN Large Hadron Collider (LHC) 
[[Bill, whose properties are found to be close to those expected for the Higgs boson, h, in the 
standard model (SM) [fJl @] [5). The rare two-photon decay of the Higgs boson has offered a 
golden detection channel in this great endeavor; for a review on Higgs phenomenology see for 
instance J6). As more events are being accumulated, its detailed properties will be uncovered. 
More optimistically, a linear collider or even a dedicated Higgs factory will be very helpful to 
reveal whether the particle is genuinely the last missing piece of SM or the first new member 
in a bigger picture. For this purpose, it is necessary to examine exhaustively its production and 
decay properties. In this work we study its multi-photon decays which are closely related to the 
two-photon channel. 

As a neutral particle, the photonic decay of the Higgs boson is a purely quantum effect. In 
SM the decay arises at the one-loop level through interactions with the vacuum charged par- 
ticles which gain mass through the Higgs mechanism (71[8l|9). This last feature implies that 
the heaviest charged particles, the top quark and the W^ 1 gauge bosons, dominate the decay 
process. The two-photon decay rate in SM was computed long time ago lH^[TTl[T2l[T3l[T4) . As 
the number of photons increases, the computation rapidly becomes technically more and more 
challenging. This arises due mainly to the constraint from electromagnetic gauge invariance, 
which results in significant cancellation amongst the contributions from individual Feynman 
graphs 031 • Fortunately, it has also been known for long |[T0l fT2| that the two-photon decay 
amplitude can be reasonably well approximated by the leading terms in the expansion of the 
Higgs mass m/, over twice masses of the relevant heavy charged particles, i.e., 2m ti w\ ion in- 
stance, for yyih = 125 GeV as indicated by the measurements at LHC, the leading terms account 
for about 65% of the exact one-loop decay width. With more photons in the final state sharing 
the Higgs boson mass and thus becoming softer on average, we expect that the expansion in 
r tt w = [ m h/(^ m t,w)] works even better. In this work we report the results for the four- and 
six-photon decay rates in this approximation. 

The low-energy effective Lagrangian for the system of the Higgs boson and electromagnetic 
fields starts with the terms that involve the least number of derivatives for a given number of 
fields while respecting electromagnetic gauge invariance. Since the Higgs field h is neutral, 
the electromagnetic field must form into gauge invariants in itself. There are two such basic 
invariants involving the least number of derivatives, & = \F^ V F^ V and Sf = ^e^ va ^ F^ v F a ^, 
which are respectively a Lorentz scalar and pseudoscalar. Attaching derivatives to h, & or £f 
does not harm the gauge symmetry but will yield terms of a higher order in the low-energy 
expansion. We assume that charge conjugation (C) and parity (P) are good symmetries, which 
is the case in SM at the level under consideration. Then, the Higgs boson as a C- and P-even 
particle can only decay into an even number of photons, and the pseudoscalar invariant must 
appear as an even power in interactions with the Higgs boson. The leading terms in the single-/? 
sector can therefore be parameterized as follows, 
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The effective coupling £y accompanies the form & l< £ } ', and is proportional to e 2 ( !+ ^y/jM 3_4 ( !+J ), 
where e is the electromagnetic coupling, is the Higgs coupling to the relevant heavy charged 
fields (top quark and W ± bosons in SM) of mass M, and a sum over them is implied. 

To present the decay amplitudes in a compact form, we introduce the notations for a pair of 
photons, ij, with outgoing momenta ki, kj and polarization vectors £ ; , £/, 

Eij = ki • £j k j • £j kj ■ kj£j • £j , (2) 

Oij = e^pk^e^ef (3) 

A factor of & (&) in Jf h when assigned to the photon pair ij will yield in the amplitude a factor 
of Ejj (0^). Note that both notations are symmetric under the interchange of the two photons. 
The amplitude for the two-photon decay h—t 12 is 
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£l £l2 (4) 



The results for the decay into more photon pairs are similar with the only difference being in 
combinatorics. For instance, for the four-photon decay h — > 1234, 

M234 = E,20E 1234 + ^02^12345 (5) 

where 

£l234 = E12E34+ E13E24 + Ei4E23, (6) 

01234 = 012034 + 013024 + 01 4 O23, ^ 

and for the six-photon decay h — > 123456, 

M2-6 = ^30^12-6 + ^12 [^1203456 +14 perms.], (8) 

where 

£12-6 = £12(^34^56 + ^35^46+^36^45) +4 perms. (9) 

The sum over photons' polarizations in the amplitudes squared is complicated, and has been 
done with the help of a computer code. We find for the two-photon decay, 
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£i 2 2£/(12), (10) 



where 

U(12) = (12) 2 (11) 

with (ij) = ki ■ kj. The pure-i? and pure-0 terms for the four-photon decay are found to con- 
tribute the same, so that 

£Ki234| 2 = [fe 2 o + S [I0t/(1234)-8X(1234)] 

+2&o&2[-6t/(1234) + 8X(1234)], (12) 
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where 



£7(1234) = £ (ij) 2 (mn) 2 = (12) 2 (34) 2 + 2 perms., (13) 

3 perms. 

X(1234) = £ (ij)(jm)(mn)(ni) = (12) (23) (34) (41) +2 perms. (14) 

3 perms. 

The result for the six-photon case is fairly lengthy: 

= £2, [44*7(1 ■ • -6) +24V(1 • ■ -6) + 120y(l • ■ -6) -56X(1 ■ • -6)] 
+^ 2 2 [348f/(l---6)-8\/(l---6)+504y(l---6)-216X(l---6)] 
+2^12 [-361/(1 • • -6) - 104V(1 ■ • - 6) -360y(l • ■ -6) + 168X(1 • • -6)] , (15) 



where the functions of six-photons' momenta are 



(16) 
(17) 
(18) 
(19) 



The phase-space integral for an all-massless final state can be completely finished for a 
polynomial of momenta. An elegant method was provided in [16]. Denoting the n-photon 
phase-space integral as 
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^'hfnh^^ir-p'l <20) 

where p is the incoming momentum of the Higgs boson, the basic idea is to decouple the 
correlation amongst the photons' momenta by converting the delta function back to an integral 
over the coordinate x and introducing a separate integral for each photon lfT6ll . 



where an infinitesimal, negative imaginary part is attached to jco to make the integral well- 
defined. The factors of k are induced by derivatives with respect to x, for instance, for our study 
here, the following will be required, 

J a P(x) = J £Ak a kPe- ik * = i 2 d a dh(x), (22) 
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and the scalar products of different photons' momenta in the phase- space integral are obtained 
by appropriate contractions with the signature tensor. The final central integral is the following 
one for m > 2 lfT6l : 

We list the integrals that are required for the Higgs decays up to six photons: 

[2r], = PS 2( 12 ) 2 = ^6 12/ 4(P), (24) 

[4y] 1 = PS 4 (12) 2 (34) 2 = ^^144/ 8 b), (25) 

(4k) 4 

[4y} 2 = PS 4 (12)(23)(34)(41) = ^i^ n 84/ 8 (p), (26) 

[6/], = PS 6 (12) 2 (34) 2 (56) 2 = ^^ 1728/ 12 (p), (27) 

(4k) 6 

[6y] 2 = PS 6 (12)(23)(34)(45)(56)(61) = ^-^ 7327,2^), (28) 

(4k) 6 

[6y] 3 = PS 6 (12) 2 (34)(45)(56)(63) = ^pl008/ 12 (^), (29) 

[6y] 4 = PS 6 (12)(23)(31)(45)(56)(64) = ^)T8 576/ 12 (p) (30) 

Noting that the terms obtained by permutation of the photons' momenta contribute the same in 
the phase-space integrals, the decay widths are worked out to be 

r(2y) = ll^o, (3D 

m 11 

r (47) = 220 33 h 52ln5 [2(& + &)-^02l (32) 
m 19 

r (67) = 2353853% ll7r 9 ^30 + 9^2-3^2] (33) 

The above results apply to any theory that induces effective interactions between a single 
Higgs scalar and various numbers of photons that preserve C and P symmetries as shown in eq. 
(OQ). In the following, we will show the numerical results in SM. To this end, we have to work out 
first the £ parameters. Since the fermions and weak gauge bosons in SM gain mass completely 
through interactions with the Higgs field which develops a vacuum expectation value, v, their 
mass terms and interactions with the Higgs boson appear in the combined form m(l +h/v) to a 
linear or quadratic power. This implies a low-energy theorem ffT2| between a process involving 
a Higgs particle and one without. For a system of constant electromagnetic field strength and 
Higgs field, we can thus establish a relation between the effective Lagrangian Jz?/, in eq. ® and 
the Heisenberg-Euler effective Lagrangian J2?h.-e.: 

*k = -( Em/^-+E^^-W"E. = — ( L m fA +E4 A W-e. (3 4 ) 

v \y omf y amy J mw \ / " m f v " m v J 
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The mass -proportionality means that only the top quark and W ± bosons are actually relevant; 
namely, for photonic Higgs decays, we only have to include in J2?h.-e. the polarization effects 
due to vacuum top quarks and W ± bosons: 



c/> cat , c/>W | 

~£H.-E. — -%-E. + -^H.-E. + 



(35) 



The vacuum polarization effects due to top quarks Jzf^ E can be readily obtained from the 
QED result [fT71fT8ll . see [fT9l for a review, by the substitutions of the electron parameters with 
the top ones, m e — > m t and —e^-Q t e with Q t =2/3, and attaching a color factor N c = 3: 



Nr 
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ds 



Q 2 e 2 F_F+ 



t<m(Q t eFs) tanh(Q t eF+s) ' 



where 



(36) 



(37) 



The vacuum polarization effects due to a charged vector field were first studied in ||20l for a 
general electromagnetic theory of a massive charged vector field. It was found that the problem 
was exactly solvable only when the charged vector has a gyromagnetic ratio equal to 2 whence 
the charge renormalization constant was found to be larger than one in contrast to QED. In 
retrospect, the first feature means that a consistent theory of charged vector fields has to be a 
non-Abelian gauge theory while the second implies asymptotic freedom in such a theory. The 
issue was later reexamined in a spontaneously broken gauge theory in unitary gauge 11271 and 
nonlinear gauge [22], see [1231 for a review on the topic and related ones. The effective 
Lagrangian was found to be 



1 r ds _ m i s e 2 F-F+ [I - 2cos(2eF_s) - 2cosh(2<?F+s)] 



167T 2 



sm(eFs) sinh(eF + s) 



(38) 



Upon expansion in the electromagnetic field, the low-energy theorem in eq. (1341) then yields 
the dominant contributions in SM to the couplings in eq. (UJ): 
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(39) 
(40) 
(41) 
(42) 
(43) 



Two features are noteworthy. First, it seems that the W ± contribution to the effective couplings 
is always negative while the top-quark contribution is alternate in sign as the power of the 
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electromagnetic invariants £f increases. This is in contrast with the impression that one 
might have got from the destructive interference between the fermionic and bosonic loops in 
the two-photon decay channel. Second, because of larger mass and smaller charge the top- 
quark contribution becomes more and more suppressed as the power of invariants increases. 
For instance, its relative contribution in the cubic terms already drops down to 10 4 and can be 
totally ignored. At m/, = 125 GeV, these parameters yield the following approximate photonic 
decay widths of the Higgs boson in SM: 

r(2y) = 5.9xl(T 6 GeV (44) 
r(4y) = 3.1x10" 15 GeV (45) 
r(6y) = 8.5xl0~ 27 GeV (46) 

We have studied the multi-photon decays of the Higgs boson that are closely related to the 
two-photon decay employed as a golden detection channel in the recent experimental discovery. 
We have parameterized the electromagnetic interactions of the Higgs boson by an effective 
Lagrangian, and relate it by a low-energy theorem in standard model to the Heisenberg-Euler 
effective Lagrangian of electromagnetism that is induced by vacuum charged particles. For the 
two-photon decay width, we recover the known result for the leading contribution that accounts 
for about 65% of the exact one-loop result. The four- and six-photon decays are found to be 
very much suppressed by phase space, in addition to naive counts of the fine- structure constant. 
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